Graphical notation reveals topological stability criteria for collective dynamics in 
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We propose a graphical notation by which certain spectral properties of complex systems can be 
rewritten concisely and interpreted topologically. Applying this notation to analyze the stability of 
a class of networks of coupled dynamical units, we reveal stability criteria on all scales. In particular, 
we show that in systems such as the Kuramoto model the Coates graph of the Jacobian matrix must 
contain a spanning tree of positive elements for the system to be locally stable. 



Discovering how the interactions between constituents 
of a system determine its macroscopic behavior is a cen- 
tral aim of physics. Ideally, we find properties in the 
system's detailed organization that have direct system- 
level implications. Significant progress has been made in 
identifying such properties on a local level. For instance 
the implications of the degree distribution, the distribu- 
tion of a network's links among the nodes, is well under- 
stood By contrast, identifying meso-scale properties 
that have a distinct effect on macroscopic behavior, de- 
spite recent advances 01, remains a difficult challenge. A 
starting point is often to link the macro-behavior to spec- 
tral properties of certain matrices, such as the adjacency, 
the different graph Laplacians, or the Jacobian matrix of 
a dynamical system. Here again, the dependence on lo- 
cal properties, i.e., certain rows or diagonal elements, can 
be discovered with relative ease. By contrast, analytical 
insights on how meso-scale patterns affect the spectrum 
are difficult to obtain, as they typically require the eval- 
uation of determinants leading to complex expressions. 

In this Letter we propose a graphical notation, rem- 
iniscent of Feynman diagrams, that facilitates comput- 
ing spectral implications of meso- and macro-scale struc- 
tures. We illustrate the usage of the notation by in- 
vestigating dynamical stability of stationary and phase- 
locked solutions in a class of symmetrically coupled dy- 
namical systems, containing for instance the Kuramoto 
model Here, the proposed notation greatly reduces 

the complexity of the mathematical expressions and en- 
ables the derivation of topological stability criteria on all 
scales. In particular, we show that in the systems under 
consideration the graph defined by the off-diagonal ele- 
ments of the Jacobian must contain a spanning tree of 
positive links to admit stable solutions. 

Any given matrix J can be interpreted as defining the 
connectivity of an abstract weighted network Q, the so- 
called Coates graph. In this network two nodes i,j are 
connected if Jij ^ 0. We can interpret a product of 
different matrix elements as defining a (not necessarily 
connected) network motif, e.g. JijJjkJki corresponds to 
a triangle path. The advantage of this graphical read- 
ing is that complex structures appearing in systems of 
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FIG. 1. Examples for the graphical notation. Symbols de- 
note the sum over all non-equivalent possibilities to build the 
depicted subgraph with the q nodes € S. Plotted are two ex- 
ample terms and their algebraic and topological equivalents. 



equations can be described using network terminology. 

It is useful to define a basis of symbols, 
x,|,A,n,0,..., denoting sums over all non-equivalent 
cycles of length n = 1,2,3,4,5,..., respectively. We 
allow concatenation of these symbols by the '■' sign. 
For a given set S of node indices, products of symbols 
denote the sum over all non-equivalent possibilities to 
realize the depicted motif with the nodes in S. For 
instance, if 5' = {«, j, fc, m} then x • □ represents the 
sum over all products of elements that would be read 
as a cycle containing four nodes and a self-loop, e.g. 
JiiJjkJkiJiviJmj (cf. Fig. [T|). Finally, we define $s as 
the sum over all products corresponding to all acyclic 
graphs that can be drawn by placing \S\ (undirected) 
links such that each link starts at a different node in S. 

In the remainder of this Letter we illustrate the usage 
of the proposed notation by considering the synchroniza- 
tion of dynamical units. This subject was chosen as it 
is presently of broad interest in physics and appears in 
many fields including biology, ecology, and engineering 
5-0I- The paradigmatic model proposed by Kuramoto 
3| opened the field for detailed studies of the interplay 
between the structure of the interaction network and col- 
lective phenomena [8-ll|. These revealed the influence 
of various topological measures, such as the clustering 
coefficient, the diameter, and the degree or weight distri- 
bution, on the propensity to synchronize 12l-ll4|. How- 
ever, recent results indicate that beside global 
topological measures also details of the exact local config- 
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uration can crucially affect synchronization. This high- 
lights synchronization of phase oscillators as a promising 
example in which it may be possible to understand the in- 
terplay between local, global, and mesoscale constraints 
on stability, that severely limit the operation of complex 
technical and institutional systems [13 • 

To provide a specific example we focus on the Ku- 
ramoto model Q , describing a network of N bidirection- 
ally coupled phase oscillators 



Xi = CJj + 



sin{xj — Xi) , Vi e 1 . . . iV 



(1) 



where Xi and Wi arc the phase and intrinsic frequency of 
oscillator i and A is the weight matrix of an undirected, 
network. The model can exhibit phase-locked states, 
which in a co-rotating frame correspond to steady states 
of HI) . The local stability of such states is determined by 
the eigenvalues of the Jacobian matrix J G R^^^ defined 
by Jik = dxi/dxk- If all eigenvalues of J are negative, 
then the state under consideration is locally stable. 

For deriving topological stability criteria we apply Ja- 
cobi's signature criterion (JSC), also called Sylvester cri- 
terion. The JSC states that a hermitian matrix J with 
rank r is negative definite iff all principal minors of order 
q <r have the sign of (—1)'^ [19[, i.e., iff 



sgn 



(^|S|,5) -(-1)'^' V5c{l,...,iV}, \S\<r (2) 



where -D151 5 := det [Jik), i,k € S. Note that the whole 
set of equations ([2|) constitutes a sufficient condition for 
stability, whereas each equation already constitutes a 
necessary condition. 

Stability analysis by means of JSC is used in control 
theory [l9| and has been applied to problems of different 
fields from fluid- and thermodynamics to offshore engi- 
neering and social networks (2ol - [23l | . An interesting prop- 
erty of the JSC is that it provides necessary conditions 
on all scales. However, application of the criterion were 
previously limited mostly to systems with few degrees of 
freedom because of the difficulties associated with com- 
puting large determinants. 

In the following we rewrite the JSC conditions in the 
notation proposed above. By direct comparison we find 
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where symbols with n > 2 appear with a factor of 2 that 
reflects the two possible orientations in which the corre- 
sponding subgraphs can be paced out. Symbols with an 
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FIG. 2. Symbolic calculation of a determinant using the 
graphical notation. Consider the minor -D4,s of the graph Q 
sketched above. If S is chosen to be the set of nodes plotted 
in grey, the terms of the fourth JSC, -D4,s, can be written as 

X-X-X-X = J11J22J33J44 = (-l)*(Jl2-|-Jl3)(J21-|-J23)(J31 + 
J32+J34)iJ43 + J45) = A + B + C + 2D, -x-x-\ = -{B + 2C), 
I ■ I = C, 2 X -A = —2D and — 2n = 0. We can immediately 
read off that D4,s = A = $s, defined in the text. 



even (odd) number of links carry a negative (positive) 
sign related to the sign of the respective index permu- 
tation in the Leibniz formula for determinants [24|. We 
note that for instance that De.s ba-s only 11 terms in the 
proposed graphical notation in contrast to 720 terms in 
conventional notation. 

The highly aggregated but somewhat awkward Eq. ([3]) 
applies to all systems of symmetrically coupled one- 
dimensional units. A more convenient expression can 
be obtained by focusing specifically on systems having 
a zero row-sum (ZRS), such that Ju = ^'^j^i Jij- In 
the Kuramoto model (and many other systems), the ZRS 
results from a force balance along the links. In the topo- 
logical reading, it allows replacing all self-loops Ju by the 
negative sum over all links connected to respective node 
i. For the first term, x''^', of a minor D^g^ g this leads 
to a sum over all graphs that can be drawn by placing 
1 5*1 (undirected) links such that every link starts from a 
distinct node in S. By elementary combinatorics it can 
be verified that all other terms of D^g^ g cancel exactly 
those subgraphs in xl'^l that contain cycles (cf. Fig. 2), 
leaving exactly <I>s defined above. For systems with ZRS, 
we can thus express the minors as 



(4) 



/|5|,5 = (-l)l'^l$s- 

Thus the JSC stability conditions translate to 

$s > 0, y S with \S\<r . (5) 
We remark that $5 with \S\ = iV — 1 is the sum over 
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all spanning trees of the network, such that Kirchhoff's 
Theorem [25|, |26[ appears as the special case of Eq. Q . 



The graphical stability conditions $s > conveniently 
conceals the complexity of the underlying determinants. 
Below we show that results can now be obtained by rea- 
soning on the graphical level, i.e., without digging up the 
complex underlying expressions. 

Consider a network containing a tree-like branch that 
is only connected to the rest of the network by a single 
link. We choose 5" as the set of nodes that are located in 
the branch and focus on the condition $5 > 0. The rules 
for constructing the network motifs in $5 imply they 
must contain one link starting from every node in S. By 
starting from the nodes of degree 1 (having only one link) 
and working downward one can see that there is only one 
motif contributing to 4>5, which consists of all links in 
the branch and the link connecting it to the rest of the 
network. The condition $5 > therefore implies that the 
number of links in $5 associated with a negative elements 
must be even. Further, if this number were greater than 
zero, one could always find a part of the branch to which 
the same conditions apply, but which contains only one 
of the negative links, such that a stability condition on a 
smaller scale is violated. Thus, stability requires that all 
links appearing in such tree-like branches correspond to 
positive entries of J. 

Using similar arguments as above the implications of 
different meso-scale motifs can be determined. The anal- 
ysis reveals restrictions on (a) the number and position 
of potential negative links and (b) the absolute value of 
their weights. We find that restrictions of type (a) can be 
subsumed under one general condition: To admit stable 
solutions, a dynamical system with symmetric Jacobian 
and ZRS must possess a spanning tree made up entirely 
of positive elements. 

To prove the statement above consider that in a net- 
work without a positive spanning tree it must be possible 
to partition the nodes into two nonempty sets Ii , I2 such 
that 



Jij < y i,j \ i e h, j e I2 



(6) 



The idea is now to evaluate the stability conditions (HI for 
different S C Ii thereby exploiting that all links leading 
out of Ii have negative weights. It is convenient to define 
E* as the set of links connecting Ii and I2, and X = 
{xi, . . . , Xm} as the subset of nodes G Ii incident to at 
least one link from E* ('boundary nodes'). Further, we 
define ai as the sum over all elements of E* incident to 
Xi, and, for any subset Y of X, ay ■— YimeY '^■m- Finally, 
for any subset y of X , we define Ty as the sum over all 
forests of Q that (i) span /i, and (ii) consists of |y| trees 
each of which contains exactly one element from Y . We 
can now write 



^ii\c = X] '^btbuc 

BCX\C 



(7) 




FIG. 3. Decomposition of a graph Q in acyclic parts (black), 
unbranched segments of cycles (grey) and branching points 
(open black). Lines represent paths of Q that may contain 
several nodes. Stability requires that (i) the acyclic parts 
only contain links with positive weights; (ii) any unbranched 
segment of a cycle contains at most one hnk with negative 
weight; (iii) the number of links with negative weight does 
not exceed the number of independent cycles (here, 5). 



where B and C are disjoint subsets of X. We show that 
this is incompatible with the stability condition ([5|) by 
the contradiction 
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Therefore the existence of a spanning tree of positive el- 
ements is a necessary condition for stability. 

Note that the (global) spanning tree criterion has dis- 
tinct implications for meso-scale properties of Q. Any un- 
branched path can maximally contain one negative link, 
and the number of negative links in the network is limited 
by the number of independent cycles (Fig 3). 

Let us now turn to stability conditions of type (b), 
which restrict the absolute value of potential negative 
links and result from smaller scale stability conditions: 
Consider for instance an unbranched segment of a cycle 
of Q that consists of d links c, with weights Wi, one of 
which, say Wx, is negative. Evaluating Eq. ^ for a series 
of sets 5* C {1, . . . ,d}, we find that stability requires 



< 



n 



iei' 



all distinct products of (d — 2) factors Wi, i £ I* ' 



where I* = {1, . . . ,d}\x. 

Applied to the Kuramoto model, the criterion derived 
above implies that in any phase locked state a spanning 
tree must exist on which the phase difference between any 
two coupled oscillators obeys |a;j— x^j < tt/2. In networks 
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where this condition is violated the criterion points to 
local interventions that enhance synchronizability. 

We note that the results derived above are not contin- 
gent on the specific form of the Kuramoto model. They 
apply to all symmetrically coupled systems obeying the 
ZRS condition, which includes all systems of the form 



Xi — Ci + 



Oij{xj - Xi) , Vi G 1 . 



.N 



(8) 



where the Aij are the weights of a symmetric coupling 
matrix and Oij are odd functions. Besides systems of 
coupled oscillators, Eq. ([8|) can for instance describe a 
range of variants of the Deffuant model of social opinion 
formation and ecological meta-population models. We 
emphasize that the derived results remain valid for het- 
erogeneous networks containing different Aij, Oij, and 

Although the ZRS simplifies the derivations above, 
similar calculations can be carried out for systems vi- 
olating the ZRS condition, e.g. We applied the 
graphical notation to an adaptive extension of the Ku- 
ramoto model, where the ZRS is violated and the con- 
nection strength coevolves with the dynamics on the net- 



work according to d/dt{Aij) = coa{xj 



'h- An. Such 



systems have recently received much interest in physics 
because of their role in neuroscience pTi - liij . Because of 
space constraints we postpone the detailed discussion of 
the model to a subsequent publication, but note that pro- 
ceeding similarly as above leads to a stronger spanning 
tree condition for the adaptive system: Stability requires 
in this case the existence of a spanning tree where linked 
nodes obey \xj — Xi\ < ir/A. 

The proposed notation is also applicable to certain 
questions not pertaining to dynamics, for instance to the 
question of isospectrality of hermitian matrices [s^, [111 ■ 
The key idea is that the characteristic polynomial x 
of a hermitian matrix A S ((^"xn expressed as 

x(A) = Dn{A — AI). The structure of the graph Q as- 
sociated to A — AI reveals the symbols contributing to 
X- One can then determine which changes of off-diagonal 
entries leave all contributing symbols and thus the spec- 
trum invariant. 

In the present Letter, we proposed a graphical nota- 
tion that facilitates the computation of spectral proper- 
ties of complex systems. Applying this notation to sys- 
tems of symmetrically-coupled one-dimensional dynami- 
cal units, we showed that any system obeying a simple 
force balance condition (ZRS) has to obey a global sta- 
bility condition: Local dynamical stability of steady (e.g. 
phase-locked) states requires that the Coates graph of the 
Jacobian has a spanning tree of positive links. This crite- 
rion is complementary to results obtained by other meth- 
ods (master stability function, ensemble simulations) and 
pertains to a large class of systems studied in physics, 
containing the Kuramoto model. Along with similar rules 



that can be derived analogously, the spanning-tree crite- 
rion has distinct meso-scale consequences, limiting for in- 
stance the number, position, and strength of negative ele- 
ments in network motifs. Beyond the coupled oscillators, 
the proposed approach is applicable to questions ranging 
from adaptive networks in neuroscience to isospectrality 
problems in condensed matter. 

We are grateful to Jeremias Epperlein and Stefan Sieg- 
mund for valuable discussions. 
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